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Abstract. We investigate equal spheres packings generated from several experiments and from a large 
number of different numerical simulations. The structural organization of these disordered packings is 
studied in terms of the network of common neighbours. This geometrical analysis reveals sharp changes 
in the network's clustering occurring at the packing fractions (fraction of volume occupied by the spheres 
respect to the total volume, p) corresponding to the so called Random Loose Packing limit (RLP, p ~ 0.555) 
and Random Close Packing limit (RCP, p ~ 0.645). At these packing fractions we also observe abrupt 
changes in the fluctuations of the portion of free volume around each sphere. We analyze such fluctuations 
by means of a statistical mechanics approach and we show that these anomalies are associated to sharp 
variations in a generalized thermodynamical variable which is the analogous for these a-thermal systems 
to the specific heat in thermal systems. 

PACS. 45.70.-n Granular Systems 45.70.Cc Static sandpiles; Granular Compaction 81.05.Rm Porous 
materials; granular materials 



1 Introduction 

Since the earliest studies of granular materials it has been 
evident that one of the key quantities which affects the 
system's properties is the packing fraction (fraction of the 
total volume occupied by the grains). It is well known since 
ancient times that different actions and different tunings of 
a given action can generate packings with different packing 
fractions. Typically, in experiments, such 'actions' consist 
in tapping the system with vertical vibrations or by shear- 
ing or by rotating the container or by pouring the grain 
in a container. In times when grain was sold by volume, 
the preparation protocol to achieve a dense packing was 
very important [1 j. This is even reported in the gospel as 
an example of good measure: "Give, and it shall be given 
to you. Good measure, pressed down, shaken together, 
running over, will be put into your lap." (Luke 6:38). In- 
deed, depending on the system handling, one can have 
variations in packing fractions up to 15% within the two 
limits p ~ 0.555 and p ~ 0.645 which are commonly ref- 
ereed as Random Loose Packig (RLP) and Random Close 
Packing (RCP) limits. For instance, by using the fluidized 
bed technique [2|.]j one can obtain packing fractions in 
the whole spectrum from 0.555 to 0.645 by varying the 
intensity of the flow pulses. 

One of the scientists who first investigated the micro- 
scopic nature of granular packing was J.D. Bernal that in 
a stream of papers concerning the "structure of liquids" 
reported some of the most important features of the struc- 
tural organization of disordered sphere packings |l|4|5j . It 



was Bernal who pointed out that disordered packing of 
equal spheres cannot overcome the RCP limit. Fascinated 
by the simultaneous simplicity and complexity of these 
systems, he asked the following question: "Science is mea- 
surement, but what is a good measure?". Indeed, he was 
aware that depending on the kind of external actions the 
system will result in different packing fractions. However, 
he also observed that for a given external driving the sys- 
tem produces configurations with very similar and repro- 
ducible packing fractions which fluctuate in a very narrow 
range of 0.5%. 

In this paper we show that the study of the fluctua- 
tions of such reproducible packing fractions can shed light 
on the origin of the RLP and RCP limits. Indeed, from 
a statistical mechanics perspective such fluctuations are a 
measure of the way in which the system is exploring the 
accessible phase-space under a given external driving. The 
study of these fluctuations gives therefore insights about 
the accessible phase-space under given constraints. Gran- 
ular materials are particle systems in which the sizes of 
the constituents are large enough such that they are not 
subject to thermal motion fluctuations. Therefore, a direct 
application of a thermodynamical theory is not straight- 
forward. However, in recent years several extensions of 
classical statistical mechanics approaches have been pro- 
posed for these systems |6I7I8I9I10I11I12I13I14I15"] . In this 
paper we use a statistical mechanics approach to relate the 
packing fraction fluctuations with changes in the system's 
structural organisation and to understand the nature of 
the structural transitions occurring at RLP and RCP lim- 
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its. In his "Bakerian lecture", Bernal explained that in 
such systems there are two fundamental questions to be 
addressed: 1) "What is the structure?" and 2) "Why has 
it got this structure?". And he resolved that the answer 
must be searched by two means: 1) Statistical Geometry 
and 2) Statistical Mechanics. By following his footsteps, in 
this paper we use a geometrical analysis and a statistical 
mechanics approach to understand what is happening at 
the two RLP and RCP limits. 



2 Materials and Methods 
2.1 Experiments 

The experimental results reported in this paper concern 
experiments from the AAS database of disorder packings 
[16J. Specifically, we investigate the six samples 'A-F' de- 
scribed in details in |17|18|19] . They are dry packings of 
acrylic mono-sized spherical beads prepared with different 
methods in a range of packing fractions between 0.58 and 
0.64. The samples B, D-F contain approximately 35000 
beads of diameter 1.59 mm. Whereas, the samples A, C 
contain approximately 120000 beads of diameter 1.00 mm. 
Polydispersities are within 0.5% and they are placed inside 
a cylindrical container with an inner diameter of 55 mm 
and filled to a height of ~ 75 mm |17|18|19j . We also re- 
port data from 12 experiments concerning glass beads in 
water prepared at packing fractions between 0.56 and 0.60 
by means of a fluidized bed technique [2 3 . Each sample 
consists of about 145000 beads of diameter 250 ± 13 jam 
placed in a cylindrical glass container with an inner diam- 
eter of 12.7 mm. 



2.2 Numerical simulations 

We generate a set of packings by means of an event-driven 
molecular dynamic simulation of hard spheres which uses 
a modified Lubachevsky-Stillinger algorithm [20 21 22 . The 
algorithm starts from random points in space and makes 
them grow uniformly into non-overallping spheres with the 
sphere positions evolving in time according to Newtonian 
dynamics. The simulation is ended when the sphere sizes 
cannot be increased any longer and a 'jammed' state with 
diverging collision rate is reached. Large expansion rates 
produce jammed configurations with low packing frac- 
tions whereas slower growth rates lead to larger packing 
fractions. With this technique, the least dense attainable 
jammed configurations have packing fractions p ~ 0.56 
which correspond to the RLP limit. On the other hand, 
for very slow rates, crystalline nuclei with large packing 
fractions (up to the limit ~0.74) can be formed. In our 
simulations, by varying the growth rate between 500 and 
0.00001, we generate jammed configurations with pack- 
ing fractions between 0.56 and 0.65. We also generate 
non-jammed configurations in the range of packing frac- 
tions between 0.1 and 0.55 by keeping the growth rate at 
0.001 and arresting the simulation once the desired pack- 
ing fraction is reached. These non-jammed systems are 



packing models of (mostly) non-touching spheres placed 
in space without overlaps. Clearly, they are not mechani- 
cally stable. They cannot be observed in experiments un- 
der gravity but they might be relevant in studies con- 
cerning micro-gravity experiments or in colloid suspen- 
sions with matching liquid density. Some simulations with 
a modified Jodrey-Tory algorithm [23] have been also per- 
formed. The algorithm starts from a set of overlapping 
spheres with repulsive interactions. It reduces overlaps, 
until all are removed, by moving spheres and gradually 
shrinking their radii [24] . All numerical simulations use a 
cubic box with periodic boundary conditions. All the anal- 
ysis of static packing properties have been performed on 
numerical samples containing 10000 spheres. Global pack- 
ing fluctuations have been studied on systems of different 
sizes (from 200 to 10000 spheres) and repeating simula- 
tions several times (200 at least per each average packing 
fraction) . 

3 Geometrical study of a structural transition 

It has been pointed out in [25] that the 'common- neighbor 
analysis of structure', first introduced by Clarke and Jons- 
son in J26j , is a very powerful method to detect structural 
organization. Such a construction consists in considering 
couples of neighbouring spheres which stay with centers 
within a given threshold radial distance and retrieve all 
the neighbors that the two spheres have in common. In 
this paper we choose a threshold distance of 1.255 sphere 
diameter. Such a distance coincides with the one used by 
Hales to individuate 'near' spheres in his recent proof of 
the Kepler's conjecture [27] and it was also used recently 
by [24128] in the geometrical study of the RCP limit. This 
is a rather convenient distance: not too small in order to be 
little sensitive to local rearrangements, and not too large, 
in order to avoid the averaging out of local properties. 
Let us remark that the choice of such threshold is not 
critical. The properties reported in this paper (Figj2| are 
consistently observed in a range of thresholds from 1.05 
to 1.4. 

We calculate the fractions p(q) of couples with q com- 
mon neighbors respect to the total number of couples. Fig- 
ure [l] shows the results of such analysis performed on 15 
numerical simulations (modified Lubachevsky-Stillinger al- 
gorithm), 6 experiments with acrylic beads in air (A-F) 
and 12 experiments with glass beads in water (fluidized 
bed technique). Only fractions with q = 2,3,4,5 (which 
have the largest statistical weight) are shown. The fig- 
ure reveals a very good agreement between simulations 
and experiments indicating that the structural properties 
are little sensitive to the preparation method and to the 
physical characteristics of the grains. Conversely the fig- 
ure reveals a clear universal dependence of the structural 
properties on the packing fraction p with opposite trends 
for £>(2), p(3) and p(4), p(5). 

A very clear signature that something is occurring to 
the structure around the Random Loose Packing pack- 
ing fraction is revealed by the sharp changes in the be- 
haviors of p{2) and p(4) occurring between p = 0.55 and 
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In order to better visualize any structural transition, 
in Fig(2]we plot the ratio: 



Fig. 1. Fraction of common neighbors calculated from: 
1) packing models of un-jammed packings of spheres 
('+'); 2) jammed packings simulated by using a modi- 
fied Lubachevsky-Stillinger algorithm [20121122] with different 
growth rates ('x' ); 3) six experiments A-F [17118119] with 
dry acrylic beads ('o, □, *, o, A, y'); 4) twelve experiments 
with glass spheres in water prepared by means of fluidized beds 
technique 3 ('>')• 




Fig. 2. The ratio A = (p(3)p(5))/(p(2)p(4)) vs. p reveals very 
sharp changes at the packing fractions corresponding to the 
Random Loose Packing and Random Close Packing limits (the 
two vertical lines). Symbols are as in Fig[l] 



p = 0.56. Some evidences of the onset of a different regime 
can also be observed from the behavior of p(3) and p(5) 
around the Random Close Packing limit (p ~ 0.645). Such 
a change must be due to the fact that in the crystalline 
phase (which, in the simulated samples, begins to nu- 
cleate above the RCP limit) there are no configurations 
with 3 or 5 common neighbors. Indeed, any closed packed 
phase, made by stacking hexagonal layers of spheres (Bar- 
low packings [27]). can admit only 2 or 4 common neigh- 
bors. 



A = 



P(3)p(5) 



(1) 



which is a 'signature' of disordered arrangements weight- 
ing configurations which are impossible in crystalline pack- 
ings (3 and 5 common neighbors) against configurations 
that instead are common in crystals (2 and 4 common 
neighbors). The plot of A reveals two sharp transitions 
occurring respectively at the RLP and at the RCP limits. 

We have therefore acquired a first evidence that there 
are neat structural changes occurring at the RLP and RCP 
limits. Now let us establish if these geometrical changes 
are associated with other changes in the statistical me- 
chanics properties of these systems. 



4 A statistical mechanics study 

As discussed in the introduction, different preparation pro- 
cedures and different experiments can result in granu- 
lar packings with different total occupied volumes V (or 
equivalently different packing fractions p = irNd 3 /(6V) 
with N the number of spheres and d their diameters) . 
Here we are interested in the properties of the set of all 
possible total volumes (packing fractions) which can be 
reached by means of a chosen system's driving. 

Any statistical mechanics theory will typically yield 
to an expression for the probability distribution of the 
volume fluctuations at equilibrium of the following form 
[3129] : 

rfl^^S^ ; (2) 



Ev W)e 



where Q(V) is the number of microscopic states which are 
classifiable under the same (coarse grained) state with vo- 
lume V. The quantity x _1 is a temperature-like intensive 
variable which was named 'compactivity' by Edwards [6]. 
It is determined by the constraint on the average volume: 



V = {V) = Y^VVoo{V) 



(3) 



A derivation of Eqj2] from a minimal set of statistical ar- 
guments is provided in [3 ; whereas a complete deductive 
statistical mechanics derivation is given in [29] . 

The challenge is to compute the number of equilibrium 
configurations £2(V) associated with states which occupy 
a total volume V under a given system preparation. To 
this end we can image that the whole system is made of 
a number k of 'elementary cells' {ci,...,Cfe} [29 3. Let 
us stress that the number of such elementary cells does 
not coincide in general with the number of grains in the 
system. Given such cellular partition, f2(V) can be com- 
puted exactly, under the two following assumptions: (1) 
these cells can have arbitrary volumes above a minimum 
value v m i ni under the sole condition that the whole system 
must occupy a total volume V; (2) all the cell-properties 
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Ci are either completely determined by their volumes Vi 
or they are independent from Vi. In this case, we have 



Q{V) 



1 \ k fV rV rV 



dv k 



5(v 1 +v 2 + ... + v k -V) 



(V - kv min ) 

A 3k (k-1)\ 



k-1 



(4) 



with A a constant analogous to the Debye length. Substi- 
tuting into Eqj2j and by using Eq(3]we obtain 

(V - kv min ) 



and 



Pco(V) = f(V,k) = 
r(k) {V-V m in) k 



exp 



V-Vmi, 
'V-Vmir, 



(5) 



, (6) 



with Vmin — kv m i n . The function f(V,k) is the probabil- 
ity density function to find a packing of k elementary cells 
occupying a volume V when the system is subject to an ex- 
ternal driving that produces an average occupied volume 
V. Note that Eqjijis valid for any k. Indeed, the observ- 
able system can be any arbitrary sub-set of a larger sys- 
tem. Moreover, the experiment can be performed either on 
several different independent systems or - equivalently - 
on several non-iteracting sub-sets of a large system. Eqj6]is 
a Gamma distribution in the variable V — V m i n ; it is char- 
acterized by a 'shape' parameter k and a 'scale' parameter 
X [30]. In Ref.[29 we named such distribution k-gamma 
distribution. For this distribution the average volume (V) 
coincides with V and the variance is 



k 



(7) 



This last relation is very useful because it provides a prac- 
tical means to evaluate k from a set of volume measure- 
ments: k = (V - V min ) 2 /a%. 



5 Granular temperature, fluctuation relation 
and Specific heat 

Equation [6] predicts that the statistical distribution of 
the volume fluctuations depends only on the parameter 
k which counts the number of elementary cells in the sys- 
tem. Let us better understand the physical and statisti- 
cal mechanics meaning of such quantity. Following Ed- 
ward's ideas [6 7 , in granular systems a 'granular tem- 
perature' (compactivity x) can be inferred from an anal- 
ogy with the thermodynamical relation (3 = 1/(&bT) = 
d(Entropy) / d (Energy) [60, by susbstituting the volume 
to the role played by the energy in thermodynamical sys- 
tems. In the present approach we can write the 'statistical 
entropy' (or Gibbs entropy) for an ergodic set Z charac- 
terized by and average volume V, as [29] : 



S(Z) 



vez 



,(V)\og Poo (V) 



5>oo0W) 

vez 



(8) 



Random 
Loose 
Packing 



Random 




Fig. 3. Behavior of the standard deviation of the Vorono'i 
volume fluctuations (a v ) vs. packing fraction (p). Very sharp 
changes are observed at p ~ 0.555 and at p ^ 0.645. Symbols 
are as in Fig[l] 



which, in the notation used in this paper, becomes 



S(Z) = k 



1 + ln 



kA d 



leading to 



A 



gr 



dS(Z) 
dV 



X 



v-v m 



(9) 



(10) 



The Edwards' compactivity x — Pgr [6 1 7] is therefore the 
average free- volume per elementary cell x = (V — V m i n )/k. 
This means that, in the present approach, the 'granular 
temperature' is a measure of the kind and the degree of 
space-partition into elementary cells. The volume fluctua- 
tions within the ergodic set can be directly calculated from 
Eqj6] and one can verify that the correct relation between 
compactivity and volume fluctuations is attained: 



x 2 ^ 1 = ((v - (V)f) 



dx 



(11) 



From this equation, substituting Eq. [7J we obtain the 
following relation for the parameter k: 



k = 



d(V) 
dx 



(12) 



The parameter k measures therefore the amount of vol- 
ume that must be added to the system in order to increase 
of one 'granular degree' the compactivity. The analogous 
quantity for molecular gasses is: dE/dT, which is the spe- 
cific heat. In analogy with ordinary thermodynamics such 
'specific heat' is expected to be sensitive to changes in the 
system's internal properties. 
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Fig. 4. Behavior of k calculated from the fluctuations of the 
Vorono'i volumes by using k = (V-Vrmnf/al (Egg) vs. pack- 
ing fraction p. Symbols are as in Fig[l] 




0.7 p 0.8 

Fig. 5. 'Granular temperature' (compactivity x) vs - packing 
fraction p. The symbols * refer to 25 numerical simulations of 
packings with 10000 spheres generated by using the Jodrey- 
Tory algorithm 23 . 



6 Changes in volume fluctuations around RLP 
and RCP limits 

We first investigate the volume fluctuations at the level of 
a single grain. For this purpose we use the Vorono'i par- 
tition where we calculate the portion of space closest to 
a grain center respect to any other centre in the pack- 
ing. In Fig. [3] we report the standard deviation a v of the 
distribution of volumes of the Vorono'i regions inside the 
various experiments and numerical samples. One can ob- 
serve that the fluctuations change abruptly in correspon- 
dence of the two RLP and RCP limits. We also observe 
that, within an overall decreasing trend, there are small 
but sizable changes at intermediate packing fractions such 
as 0.58 ad 0.6. Such behavior is reflected in the value of 
the parameter k = (V — Vmin) 2 1&1 (Eqj7]) as reported in 
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Fig. 6. Rescaled standard deviations of the global fluctuations 
of the whole sample volume <j v (N) / for various system sizes 
(N) calculated from several thousands numerical simulations 
by using the Lubachevsky-Stillinger algorithm 20 21 22]. All 
data rescale on the same values, deviations are observed only 
for the smallest sample (N = 200, +) near the RCP limit 
where the small size of the sample increases the probability 
of crystallization. The symbols refer to standard deviation 
of the local Vorono'i volume fluctuations a v (same data as in 
Fig§. 



Fig(4j Note that the minimum volume of a Vorono'i re- 
gion in packings of equal sphere with unit diameters is a 
fixed value corresponding to the volume of a dodecahedral 
region: V m i n = 0.694... [3 . The value at zero packing frac- 
tion (k = 5.586) was calculated analytically for random 
Poisson points in three dimensions j31|32j . Figure [4] shows 
that at low packing fractions, for non jammed configura- 
tions, the value of k increases almost linearly with p. Then, 
it drastically decreases to values between 11 and 15 when 
the system gets into jammed configurations. The inset in 
the figure shows that there are differences in the values of 
k for different systems and within the same system at dif- 
ferent packing fractions. One can also note a rather sharp 
change in the experimental data occurring around p ~ 0.6 
which might indicate some kind of transition at this pack- 
ing fraction. Above the packing fraction ~ 0.645 (RCP), 
the packings contain partially crystallized regions and the 
change in the kind of structural organization is reveled by 
a sharp drop in the value of k that eventually will go to 
zero at the crystalline limit (p = 0.740...). These data are 
consistent with Refs. [3 29 where we have shown that the 
volume distribution of the Vorono'i regions follows remark- 
ably well the theoretical prediction f(V,k) (Eq. |6| with k 
in the range between 9 < k < 25. This implies that, de- 
pending on the kind of system in exam and on the packing 
fraction, there are between 9 to 25 elementary cells which 
are in average contributing to the volume of each Vorono'i 
region. The impressive fact is that these systems are very 
different (ideal Newtonian spheres, acrylic beads in air 
and also glass beads in water [3]) and they are prepared 
in very different ways (pouring, tapping, fluid flows, shear- 
ing, hard-spheres molecular dynamics). The fact that all 
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these distributions follow the same law f(V, k) (k-gamma 
distribution, Eqj6| suggests that there are universal prop- 
erties that determine the packing configurations and their 
fluctuation laws. On the other hand, the fact that dif- 
ferent systems or different preparation methods yield to 
distributions with different values of k indicates that this 
quantity is an important parameter to control and char- 



acterize the system's properties. Equation 10 reveals a di- 
rect relation between the parameter k and the Edwards 
compactivity %. The variation of the compactivity with 
the packing fraction is reported in Fig(5l One can observe 
that x nas a similar behavior to a v (Fig|3| revealing sharp 
peaks within an overall decreasing trend. Again we ob- 
serve large changes occurring at the RLP and RCP limits 
indicating that strong changes in the system's properties 
are happening at these limits. In the figure there are also 
reported data for packing models generated by using the 
Jodrey-Tory algorithm [23]. With this algorithm we can 
reach larger packing fractions spanning a region above the 
RCP limit where the system becomes polycrystalline. 

We also investigate the volume fluctuations at the level 
of the whole sample. In this case we can only study jammed 
configurations above the RLP limit. In Ref.[29] we demon- 
strated that the distribution of the total volume occupied 
by the packed spheres follows accurately well the predic- 
tion of EqjBJ The parameter k calculated from the global 
fluctuations reveals a clear peak at the RCP limit (see 
Fig. 6 in [29]) wich is consistent with the abrupt changes 
at RCP observed in the local k (FigJiJ). We observe that 
the global standard deviation of the volume distribution 
in a packing of N spheres (a v (N)) scales with the system 
size accordingly with the law: cr v (N) = (j\^/N . This scal- 
ing law is clearly demonstrated in Fig JH] where a v {N) / y/N 
for various system sizes between N = 200 to N = 10000 all 
collapse onto a single trend. Such a scaling confirms that 
the compactivity, calculated from the global fluctuations 
(Eqsj7| and 10) x = a v{N)/(V — V m i n ), is indeed an in- 
tensive parameter. However, Fig (6] reveals that the scaling 
factor a\ does not coincide with the observed variance at 
the level of a single grain a v . Such a discrepancy must be 
consequence of correlations between neighboring Voronoi 
regions [33] . As consequence the compactivity measured at 
local level is different from the one measured from global 
fluctuations. In this respect, the 'proper' compactivity is 
the one associated to the global volume fluctuations; the 
local measure is an 'effective compactivity' [29]. Intrigu- 
ingly, they coincide at the RCP limit. 



7 Conclusions 

In this paper we provide two independent evidences of 
transitions occurring in sphere packings at the Random 
Loose Packing and Random Close Packing limits. 

(i) The first evidence is geometrical and it is acquired 
from the behavior of the fraction of common neigh- 
bours. In particular, we have observed that the frac- 
tion of couples with 3, 5 and 2 or 4 common neigh- 
bours (Eqjl]) has sharp changes occurring both at the 



RLP and RCP limits. Such changes indicate that re- 
arrangements towards a different packing organization 
are occurring at these limiting packing fractions, 
(ii) The second evidence is acquired from a statistical me- 
chanics study concerning the fluctuations in the local 
Voronoi volumes. We have found that all the samples 
investigated follow the same kind of statistical distri- 
bution /(V, k) (k-Gamma distribution, Eq.[6| but they 
are characterized by different values of the quantity k. 
The value of k sharply decreases when the packing frac- 
tion crosses the RLP or RCP limits. This corresponds 
to sharp freezing of some degrees of freedom associ- 
ated with changes in the packing' s organization. We 
have discussed that the quantity k is analogous to the 
specific heat in ordinary thermodynamics. 

These evidences clearly demonstrate that transitions are 
occurring at the two RLP and RCP limits. However, it 
rests unclear whether such transitions can be described 
as 'proper' phase transitions in a statistical mechanics 
framework. This difficulty is intrinsically associated to the 
fact that the system is disordered, there are no symme- 
tries to break and we cannot introduce an 'order param- 
eter' to simply describe the structural changes. Never- 
theless, we have clearly shown in this paper that abrupt 
changes in the structural properties are occurring and 
that such changes are associated with freezing of degrees 
of freedoms and consequent contractions of the available 
phase-space. Specifically, at the RLP limit the system 
undertake important changes passing from a 'compress- 
ible gas' -like behavior to a 'rigid solid' -like behavior. If 
we constraint our analysis to mechanically stable struc- 
ture (rigid or 'jammed' structures only) then we see that 
both RLP and RCP limits are associated with exhaustion 
of realizable packings. Experimental preparation meth- 
ods typically fail to find disordered packings with den- 
sities above Prcp — 0.645 even if there exists a large 
class of layered packings (Barlow packings [27]) with pack- 
ing fraction p = 0.740.... A study of the system entropy 
(Eq|9]) reveals that the number of accessible configura- 
tions decreases approaching the RCP transition but it be- 
comes of the order of one only at the estimated 'Kauz- 
mann density' px ~ 0.66 [24] which is larger than the 
observed Prcp- Similarly, in the RLP case mechanically 
stable structures are not discovered in disordered arrange- 
ments below Prlp — 0.555. On the other hand, below 
the RLP limit jammed packings can be obtained, but 
only in special crystals with self-avoiding "tunnels" [34] . 
This seems to indicate that outside the RLP-RCP lim- 
its there might be configurations but they are isolated 
regions or points in the phase-space and they cannot be 
simply reached from small improvements on the known so- 
lutions, they have therefore infinitesimal probabilities to 
be discovered. 
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